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ABSTRACT 


An  analytical  approach  is  used  to  determine  general  results  for 
the  effects  of  both  receiver  and  source  motions  on  a CW  signal 
transmitted  through  a deep  ocean  channel  at  short  ranges.  A 
bilinear  sound-speed  profile  is  used.  The  receiver  and  source 
are  restricted  to  the  surface,  and  only  SRBR  rays  are  relevant. 
Time-dependent  expressions  for  the  total-field  amplitude  and 
phase  are  developed  for  appropriately  limited  time  intervals, 
and  numerical  results  are  presented.  General  analytical  expressions 
for  the  total  field  are  derived  and  demonstrated  to  approximate 

closely  numerical  results.  These  expressions  provide  the  basis 
for  a study  of  the  acoustical  effects  of  varying  motion  parameters 
and  initial  range.  It  is  demonstrated  that  effects  of  differences 
in  range  on  total-field  phase  rate  and  the  time  interval  between 
amplitude  maxima  are  significant  at  short  ranges  and  diminish  as 
range  increases.  Effects  on  total-field  due  to  receiver  motion 
are  shown  to  be  both  significant  and  widely  varying,  depending  on 
receiver  and  source  directions  and  speeds. 
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INTRODUCTION 

J 

Until  recently,  most  studies  of  acoustic  transmissions  in 

an  underwater  environment  have  been  limited  to  sources  and 

receivers  in  fixed  positions.  Investigations  of  energy  transmitted 

from  moving  sources  usually  have  been  limited  to  the  area  of  radar 

1 2 

and,  in  acoustics,  to  classical  treatments  ' of  line-of-sight 
transmissions  from  sources  moving  in  an  isospeed  medium.  How- 
ever, the  results  are  not  generally  applicable  in  the  ocean, 
because  such  effects  as  refraction,  boundary  reflection,  and 
multipath  summation  have  not  been  included. 

In  the  past  few  years,  however,  studies  have  begun  to  address 
moving-source  problems  in  more  realistic  oceanic  environments, 
as  exemplified  by  References  3 through  10.  Most  of  these  recent 
investigations  of  relative  motion  between  receiver  and  source 
have  been  restricted  to  a receiver,  or  receivers,  fixed  in  a static 
medium.  Studies  of  the  effects  of  receiver  motion  are  not  common 
in  the  open  literature  despite  practical  uses  including  a co- 
located source  and  receiver  in  active  moving  systems  to  form 
synthetic  aperture  arrays^  and  to  conduct  echo-sounding  measure- 
ments.^ In  part,  this  may  be  due  to  the  classical  discussions 
mentioned  earlier,  which  show  that  some  acoustical  effects  of 
receiver  motion  (e.g.,  Dopplcr-shift  frequency)  can  be  deduced 
from  an  equivalent  fixed  receiver-moving  source  problem,  assuming 
line-of-sight  transmission  in  an  isospeed  medium.  For  dynamic 
and  more  complicated  media,  however,  it  is  not  correct  that  effects 
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of  receiver  motion  can  be  determined  by  extending  results  of 
fixed  receiver-moving  source  investigations. 

The  primary  purpose  of  this  paper,  in  contrast  to  those 
above,  is  to  apply  ray  theory  in  an  analytical  investigation 
of  the  effects  of  both  receiver  and  source  motion  on  the 
transmission  of  a CW  signal  in  a deep  ocean.  In  addition,  we 
will  examine  these  effects  over  rather  short  ranges,  since  to 
our  knowledge,  the  literature  pertaining  to  acoustical 
investigations  of  relative  receiver-source  motion  is  limited 
virtually  to  relatively  long  ranges.  Our  analysis  will  permit 
a specific  description  of  receiver-motion  contributions  to  the 

total  acoustic  field,  which  necessarily  include  both  temporal 

4 

and  spatial  fluctuations.  We  assume  for  simplicity  that  the 

ocean  channel  has  horizontal  boundaries  and  a sound-speed  profile 

that  is  bilinear  with  depth,  as  shown  in  Fig.  1(a).  Further,  we 

consider  the  case  where  surface  and  bottom  sound  speeds  are  equal, 

and  where  both  the  source  and  receiver  are  located  at  the  surface. 

Therefore,  only  SRBR  (surface-reflected/bottom-reflected)  rays 
12 

exist.  In  addition,  the  motions  of  both  receiver  and  source 
are  limited  to  horizontal  paths  at  constant  velocities  denoted 

■f 

by  v^  and  v^,  respectively.  Relaxation  of  the  above  restrictions 
could  be  made,  if  desired,  at  the  expense  of  relative  simplicity. 

In  Sec.  I,  we  transform  the  moving  receiver-moving  source 
problem  into  an  equivalent  problem  in  which  all  motion  is  viewed 
relative  to  a moving  reference  frame  located  on  the  receiver. 
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Expressions  for  the  geometry,  travel  time  and  spreading  loss  of 
individual  rays  are  derived,  and  equations  for  the  total  acoustic 
field  at  the  receiver  are  developed.  In  Sec.  II,  expressions  for 
total-field  amplitude  and  phase,  valid  for  appropriately  limited 
time  intervals,  are  determined  and  provide  the  basis  for  the 
subsequent  calculation  of  numerical  results.  Examples  of  such 
results  at  both  the  minimum  (2km)  and  maximum  (10  km)  of  the 
assumed  short-range  interval  are  presented  in  Sec.  III.  Then, 
approximate  analytic  expressions  for  total-field  amplitude  and 
phase  are  derived  which  are  used  in  Sec.  IV  to  explain  the 
effects  on  total-field  of  the  variation  of  motion  parameters  and 
receiver-source  range  at  some  initial  time.  Sec.  V is  a summary 
of  the  principal  findings  in  this  paper. 

I . GENERAL  DEVELOPMENT 

We  assume  a deep-ocean  bilinear  sound-speed  profile,  as 

shown  in  Fig.  1(a),  with  equal  bottom  and  surface  sound  speeds 

c^.  This  is  a reasonable  model  in  many  ocean  channels  of 

approximately  5 km  depth.  Also,  we  choose  initially  a fixed, 

left-handed  coordinate  system,  with  origin  at  the  surface, 

where  asterisks  designate  spatial  variables  in  this  fixed 

reference  frame.  The  vertical  coordinate  z*  is  measured 

positive  downward,  and  the  minimum  sound  speed  c is  at  z*  = h , 

mm 

the  depth  of  the  SOFAR  axis.  The  sound-speed  gradients  above  and 
below  this  axis  are  denoted  by  and  g^,  respectively. 
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As  suggested  in  Fig.  1(b),  both  receiver  and  source  -X 
are  on  the  surface  and  move  on  horizontal  linear  paths  at  constant 
speeds  v^  = |v^|  and  - Iv^l/  respectively.  The  symbols  O'  and 

A 

0 denote  the  positions  of  ^ and  -X  at  a time  t arbitrarily  designated 

as  zero,  while  R is  the  distance  between  O'  and  0.  For  simplicity, 
o 

y*  is  measured  positive  in  the  O'O  direction,  while  receiver  and 
source  path  directions  are  described  by  defining  the  angles  Xq 
and  yq  measured  positive  clockwise  from  the  y*  direction  to  the 
directions  of  motion  of  ^2  and  X , respectively.  At  a generic 

/v 

source  time  t,  the  moving  source  -X  emits  a CW  signal  which 
travels  in  a vertical  plane  and  arrives  at  the  moving  receiver 
at  some  subsequent  receiver  time  t.  The  positions  of  ^2  and  X at 

A A 

times  t and  t for  a particular  time  t and  for  a particular  ray 


are  shown  in  Fig.  1(b).  Consistent  with  our  notation  for 
spatial  variables,  we  will  use  asterisks  to  denote  functions  in 
the  fixed  x*,  y*,  z*  reference  frame.  Also,  functions  of 
source  time  t will  be  annotated  with  a caret,  while  functions  of 
receiver  time  t will  have  no  such  superscript. 

From  the  geometry  and  fixed  frame  of  Fig.  1(b),  we  can  write 
an  expression  for  the  horizontal  distance  traveled  by  a ray 


emitted  at  source  time  t and  arriving  at  at  time  t as  follows: 

A A A A A O A 

R*  (t)  = {R  + (v,t)  4-  (v.t)  - 2R  v,  t cosX  + 2 R v.t  cosy 

o 1 2 ol  o o2  O 


- 2 (v  t)  (v  t)  cos  (X  - y )} 

12  O O 


We  use  the  notation  R*(t),  in  lieu  of  R*(t,t),  since  t is  in  fact 
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a function  of  t.  This  point  will  be  discussed  in  detail  soon. 

Since  both  the  receiver  and  source  are  on  the  surface, 
an  SRBR  ray  will  cross  the  SOFAR  axis  2N  times,  where  N is  a 
positive  integer  corresponding  to  the  number  of  bottom  reflections. 
We  denote  the  inclination  angle  of  a ray  by  0,  where  0 is  measured 
positive  downward  from  the  horizontal  source-receiver  direction 
in  the  vertical  plane  of  the  ray,  as  shown  in  Fig.  1(b).  Also, 

A 

we  denote  the  ray  launch  angle  at ^ by  0,  (N,t),  since  0 is 

o 

A 

completely  described  by  the  values  of  N and  t and  where  the 

subscript  b corresponds  to  the  surface  sound  speed  c^. 

5 6 A ~ 

As  discussed  in  previous  papers,  ' if  T*(N,t)  represents 

A 

the  travel  time  of  the  Nth  ray  emitted  from  ^ at  time  t,  then 

A 

source  time  t and  receiver  time  t in  Eq.  (1)  are  implicitly 
related  by 


As  suggested  in  Fig.  2,  relative  to  A , appears  to  move  in 


the  direction  of  the  vector  v = V2  ~ vi  a constant  speed 

v = | v|  , and  the  medium  appears  to  have  a spatially  uniform 
"current"  which  we  denote  by  vc  - - v^.  The  path  of  -8.  is  now 
described  by  the  angle  \ I>q,  measured  positive  clockwise  from  the 
O'O  direction  to  the  direction  of  v.  The  horizontal  angle  X (t) 
describes  the  direction  of  the  current  (or  equivalently,  the  direction 
of  receiver  motior  relative  to  the  receiver-source  direction  at 

A 

time  t)  and  is  measured  positive  clockwise  from  the  current 
direction  to  the  receiver-source  direction  at  time  t.  The  direction 

A A A 

of  the  current  at  t = 0 is  denoted  by  Xq,  horizontal  range  R(t) 

A A 

and  cos  X (t)  at  any  source  time  t are  given  by  Eqs.  (3a)  and  (3b), 
respectively,  and  the  relationships  between  model  parameters  are 
given  by  Eqs.  (3c)  and  (3d) : 


A A 2 "'I  A l /2 

R ( t)  = (R  + (vt)  + 2(vt)R  cos  ij/  ] , 

o o o 

A A 

cos  X (t)  = [R  cos  X + vt  cos(iJ>  - X )] 
o o o o 


(3a) 


2 (vt)  R cos  ] 1//2  , 

o o 

(3b) 

- 2v  v cos  (X  - Y ))1//2 
12  o o 

f 

(3c) 

and 


-1 


\i;  = tan  [v_sinY  - v.  sinX  ] [v^  cosy  - v,  cosX  ] 

o 2 o 1 o 2 o 1 o 


-1 


(3d) 


where  the  multiple-valued  function  \pQ  is  restricted  to  the 

interval  0 < ij>  < 2u. 

— o 

Since  the  medium  may  now  appear  to  move  with  respect  to 
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the  reference  frame  at  (L  , it  will  be  useful  to  distinguish 
between  expressions  valid  for  moving  and  stationary  media. 

Accordingly,  we  will  identify  expressions  valid  in  a static 
medium  with  a subscript  "s"  (although  tfl  and  may  be  in 
motion) , and  expressions  for  which  receiver,  source,  and  medium 
are  motionless  with  a subscript  "o"  (i.e.,  <fl  and  are  fixed 
at  a horizontal  distance  Rq  apart  in  a still  medium) . Thus, 
expressions  without  subscripts  are  valid  in  a moving  medium. 

: 

Ref.  12  contains  rather  general  expressions  for  launch 
angle,  travel  time,  and  spreading  loss  in  a stationary  medium 
having  a bilinear  sound-speed  profile.  With  the  prior  assumptions 
in  this  paper,  Eqs.  (3a)  and  (3b)  and  the  results  of  Ref.  12  can 
be  used  to  deduce  required  time-dependent  equations  appropriate 
to  this  problem.  The  launch  angle  of  one  of  our  SRBR  rays  is 
given  by 

tan  0.  (N,t)  = N^'1  + g'1)  (c, 2 - c 2)  c'1  [R  (t)  ] _1 
b 1 2 b m b 

•l  A A 

- [4N(g1  + 92)cb]  g1g2R(t)  . (4) 

Further,  the  appropriate  time-dependent  expression  for  travel 

time  in  a stationary  medium  is 

-1  -1  -1 
T (N,t)  - 2N(g  + g ) log[(c./c  )(l+sin6  ) (l+sin0,  ) ),  (5) 

s 12  b m m b 

where  "log"  represents  the  natural  logarithm,  0^  is  the 

inclination  angle  of  a ray  at  the  SOFAR  axis,  and  sin  0.  and 

b 

sin  0^  can  be  derived  from  Eq.  (4)  and  Snell's  law,  respectively. 

• - ' J 
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Extensive  analyses  ' of  ray  geometry,  travel  time,  and 
spreading  loss  have  been  made  for  a uniformly-moving  ocean 
channel  with  horizontal  boundaries  and  a bilinear  sound-speed 
profile.  Using  these  results  and  Eq.  (5) , we  can  write  a time- 
dependent  expression,  appropriate  for  this  problem,  which  relates 
travel  time  in  a moving  channel  to  the  travel  time  in  the 
absence  of  a current  by 

A A A A A A ^ 

T (N,  t)  = T (N,  t)  [1  + M cos  0,  (N , t ) cos  A (t)  ] , (6) 

s c b 

->  -1  ~ 
where  M = |v  I c and  both  the  launch  angle  0^(N,t)  and  the 
c 1 c‘  b ^ b 

horizontal  angle  between  the  current  direction  and  the  plane  of 
a ray,  A (t) , are  time-dependent.  To  consider  spreading  loss, 
we  assume  that  a point  source,  moving  with  velocity  v and  emitting 
a unity-amplitude  omnidirectional  CW  signal,  travels  through  a 
uniform-current  medium  with  velocity  vc  in  which  an  omnidirectional 
point  receiver  is  fixed.  Neglecting  the  effects  of  attenuation 

/\  A 

and  scattering,  the  amplitude  A (N, t)  of  a ray  at  will  be 
affected  only  by  spreading  and  boundary  losses.  Reference  15 
developed  the  relationship  between  spreading  loss  in  a uniform 
current  and  spreading  loss  in  a stationary  medium,  if  the  sound 
speed  varies  linearly  with  depth.  The  results  of  Ref.  15  may  be 
applied  above  and  below  the  SOFAR  axis  of  this  problem,  and  we 
may  write 

A A A A 

L (N,  t)  = L (N , t)  [1  + 0 (M  )]  , (7a) 

b c 

where 

■ 

" • ~~ — ■ 
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" ~ "2  2 

L_(N,t)  = R (t)  tan  6 csc0csc0  (7b) 

S b b m 

is  obtair  .d  by  applying  Eq.  (3a)  and  the  assumptions  of  this 

paper  to  the  results  of  Ref.  12. 

A A 

Denoting  the  loss  at  a single  bottom  reflection  by  B(N,t), 
the  amplitude  at  the  receiver  of  the  arrival  with  N bottom 
reflections  is  taken  to  be 


A(N, t) 


(L (N , t)  ] 


-1/2 


~ ~ N 

[ B (N , t)  J , 


(8a) 


since  scattering  losses  are  being  neglected.  The  phase  of  the 

12 

Nth  arrival  can  be  expressed  as 


<|>(N,t)  = ojT  (N , t ) - (N-l)ir  - NS(N,t) 


(8b) 


where  u>  represents  circular  frequency  in  rad  sec  at  the 

! 

source,  and  frequency  shift  due  to  Doppler  is  contained  in 

A A 

the  time-dependent  phase  terms  wT  and  S.  The  second  term  on  the 
right  side  of  Eq.  (8b)  follows  from  assuming  a ir-rad  phase  shift 

/\  /v 

at  each  surface  reflection,  and  the  term  S(N,t)  represents  the 
shift  at  a bottom  reflection. 

Combining  Eqs.  (2a),  (2b),  (4) -(8)  and  dropping  the  asterisk 

notation,  since  Eqs.  (4) -(8)  are  written  with  respect  to  the 
moving  frame  at  A , we  can  write  the  Nth-ray  contribution  to 
the  total  field  at  the  receiver: 


A A A A 

V(N , t)  = A (N,  t)  sin  [wt  - (N-l)ir  - NS(N,t)] 

A A A A /\  A 

= A(N,  t)  sin{w[t-T(N,  t)  ] - (N-I)tt  - NS(N,t)} 

= A(N,F(N,t)  ] sin  (lot  - ojT  fN , F (N , t)  ] - (N-I)tt  -NS  [N , F (N , t)  ] } 


= A (N,  t)  sin  (o)t  - 4>(N,t)] 


9 


(9a) 
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where  $(N,t)  = wT(N,t)  - (N-l)ir  - NS(N,t)  is  the  phase  in  radians 
and  A(N,t)  is  the  amplitude  of  the  Nth  ray.  The  total  field  is 
formed  by  summing  the  arrivals  and,  in  terms  of  receiver  time  t, 
may  be  written  in  the  form 


CL(  t)  sin  [tot  - (t)  ] = I A (N,  t)  sin  [wt  - 4>(N,t)]  . (9b) 

N 

For  the  short  ranges  in  this  paper,  the  above  summation  begins 
at  N = 1,  and  although  the  number  of  terms  is  countably  infinite, 
rapid  convergence  due  to  bottom  reflection  loss  permits  termination 
at  some  N = N'.  Using  Eq.  (9b),  the  total-field  amplitude,  Ct(t) , 
and  phase,  $(t),  are  given  by: 


CL  (t)  = [E  A(N,t)sin  (N,t)]2  + [£  A(N,t)cos  (N,t)]2  , 
N N 

sinMt)  = CL  * I A(N,t)sin<)>(N,t)  , 

N 

and 


(10) 

(Ha) 


cos4>  (t)  = CL  1Z  A (N , t)  cos<HN,t)  . (lib) 

N 

Because  of  the  time  dependence  of  4>(N,t)  in  the  argument  of  the 
sine  terms,  Eq.  (9b)  represents  a signal  with  slight  frequency 
spreading,  and  is  not  exactly  single  frequency.  However,  we 
elect  to  interpret  the  total  field  as  a CW  signal  of  frequency  w, 
with  time-dependent  phase  $(t),  as  well  as  amplitude  CL(t) . 


II.  TIME-DEPENDENT  APPROXIMATIONS 

The  total-field  expressions  at  the  end  of  Sec.  I are  formal 
at  this  point,  because  CL  and  ♦ are  not  known  as  explicit  functions 
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of  receiver  time  t.  Fortunately,  this  difficulty  may  be  overcome 
by  considering  sufficiently  small  time  intervals.  In  the  ensuing 
approximate  development,  we  wish  to  retain  only  those  effects 

of  source  and  receiver  motion  which  result  in  corrections  larger 

t 

than  0 (M)  to  total-field  amplitude  and  phase,  where  M = vc^  . 

Therefore,  it  is  necessary  to  retain  only  corrections  larger  than 
0 (M)  to  the  amplitude  A(N,t)  and  phase  <Ji(N,t)  of  a single  ray 
(see,  for  example.  Ref.  16)  . 

In  order  to  express  Ct  and  4 explicitly  in  terms  of  t,  the 
quantity  F(N,t)  is  required.  Since  an  explicit  expression  for 
F (N, t)  generally  cannot  be  obtained,  we  proceed  to  derive  an 
approximation  to  it  by  developing  first  a time-dependent  expansion 

A A 

for  T(N,t) , given  in  Eq.  (6) . We  begin  by  using  Eq.  (3a)  to 
expand  Eq.  (4)  to  second-degree  terms  in  source  time  t: 

tan8  (N,t)  = tan©,  - (D  R + D_R  ) (cosi|»  ) (c, /R  ) Mt 
b bo  lo  2o  obo 

+ (1/2)  (D.R  ^(3cos^i|>  -1)  + D R cos^iji  } 
l o o 2 o o 

X <cb/R0>2(Mt)2  , (12) 

where 

°x  - Ntg'1  t g/McJ  - c^c;1  , (13a) 

D2  = I4N(g"1  + 9~21)ch)~1  , (13b) 

and 

tan8.  = D R1  - D R . (13c) 

bo  1 o 2 o 

Similarly, 

~ . — — A 


and  0^o  represents  the  launch  angle  of  a ray  emitted  from  a 
source  fixed  a distance  Rq  from  a fixed  receiver.  Similarly, 

A 

Eq.  (3b)  may  be  used  to  expand  cos  A (t)  to  second-degree  terms 
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G,  = - T (N) fcos8  sin*  sin4<  + (N(g  1 + g *)  (c2  - cZ)  (c  R ) 
1 o v bo  oo  1 2 b m b o 


-1 


+ R/[(4N)(g,  1 + g_^)c  ] Jcos^B,  sin8.  cosX  cosij > ) 

o 1 2 b bo  bo  o o / 


- cos2e,  CSC0  (sine  - sin8  ) [2N  (g  *+g  *)  ] {N  (g  ^g.,1)  (c2-c2) 
bo  mo  mo  bo  12  12  bm 


+ R / [ (4N)  (g  ^+g  "Sc,]}  cosX  cosiji 

O D O O 


(19b) 


and 


G_  = cos8.  csc0  (sin0  - sin0,  ){2N(g n 1+g -1)  ) 
2 bo  mo  mo  bo  1 2 


x {N(g11+921) (cb"cm)  (cbRo)  1 + R0/[(4N)(9  1+g  1>cbJ}  ' 


(19c) 


and 


a = H(c  r“1)2M2  , 

3 bo 


(20a) 


where 


H = [J2I2  + J1I21 [2N(g11  + g21)] 


(20b) 


In  the  above. 


-1  2 

I = - (D,  R + D_R  )cot0,  cos  0,  cos\J>  , 
1 1 o 2o  bo  bo  o 


(20c) 


I2  = (3/2)  (D^1  + D2Rq)2  cos40bocos2^o  + E2cot0bocos20bo,  (20d) 


J,  = - sin0,  c sc 0 sec  6,  (sin0  - sin8_  ) , 

1 bo  mo  bo  mo  bo 


(20e) 


and 
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J_  = [-(l/2)sinV  (l-sinS  ) + (l/2)sinV  sin20  (20f) 

2 bo  mo  bo  mo 

x (1  + sin0  - 2 sin0  + 2 sin20,  ) - 2 sin0,  sin0  ] 
mo  bo  bo  bo  mo 

•3  4 

* [sin^e  cos  0,  ] 

mo  bo 

We  remark  that  the  exact  expression  for  (Eq.(20a))  contains 

2 

a term  proportional  to  M M which  we  have  dropped,  since  its 

c 

contribution  is  small  compared  to  that  of  the  other  terms  in 

A 

Eq.  (17),  assuming  that  t is  suitably  restricted.  Further, 

T (N)  = 2N(g“1+  g”1) log  fee  1 (1+sin©  ) (l+sin0^  )_1]  (21) 

o i 2 b m mo  bo 

is  the  travel  time  of  the  Nth  ray  emitted  from  a stationary 

source  at  a horizontal  distance  Rq  from  a stationary  receiver 

and  traveling  through  a motionless  medium. 

To  express  results  in  terms  of  receiver  time  t,  we  substitute 

Eq.  (17)  into  Eq.  (2a)  and  obtain 

A 2 A 

a3fc  + (a2  + 1)fc  " (fc  ~ ax)  = 0 * (22a) 

from  which 

t = F (N, t ) = (2a3)'1{-(l+a2)  + [(l+a2)2  + 4a3 ( t-c^) ) 1/2} . (22b) 

A A ^ 

We  have  demonstrated  numerically  that, for  t not  too  large,  a3t 
is  small  and  Eq.  (22a)  may  be  written  as 

A A A A 

T'  (N'fc)  + G (t)  = 0 , (22c) 

A A A 2 

where  e(t)  = a3t  represents  a small  perturbation  to  the  unperturbed 
quantity 

A A A 

(N,t)  = (a2+l)t  - (t  - a;L)  = 0 . 


(22d) 


» 
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Equation  (22d)  provides  the  unperturbed  solution  t = (t-a^)  (a2+l) 

A 

Since  e <<  1,  the  difference  between  this  result  and  the 

A* 

perturbed  solution  for  t given  by  Eq.  (22b)  must  be  small,  which 
holds  only  if  the  positive  sign  is  chosen  for  the  radical  in 


Eq.  (22b) . With  this  choice  of  sign,  if  Eq.  (22b)  is  substituted 
into  Eq.  (17)  and  the  result  is  expanded  binomially,  we  obtain 
the  following  expression  for  travel  time  in  terms  of  receiver 
time  t: 

T (N,  t)  = t + (c^  - t)(l+a2)_1  + a3(ai-t)2(l+a2)"3  . (23) 

Also,  using  Snell's  law  and  Eqs.  (7b)  and  (12),  we  have  confirmed 
that 

L_1/2  (N,t)  = L_1/2  (N) [1  + O (M) ] , (24) 

s o 

where 

L (N)  = R2  tan20,  csc0,  csc0  (25) 

o o bo  bo  mo 

is  the  spreading  loss  for  the  Nth  ray  with  the  source  and  receiver 
fixed  in  a stationary  medium. 

For  the  ranges  being  considered  in  this  paper,  it  is 
reasonable  to  assume  that  the  loss  at  each  bottom  reflection  of 


a ray  is  constant  and  that  the  corresponding  phase  shift  is 
17  18 

zero.  ' Denoting  the  constant  bottom  loss  by  Bq,  and  keeping 
only  corrections  in  amplitude  larger  than  0(M),  we  may  now  write 
the  following  expressions  for  the  amplitude  and  phase  of  the  Nth  ray: 


A (N, t)  = A (N) [1  + O (M) ] = L~1/2(N)BN 

o o 


(26a) 
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<J>(N,t)  = <x)T(N,t)  - (N-I)tt  , 


(26b) 


where  T (N, t)  and  Lq(N)  are  given  by  Eqs.  (23)  and  (25), 
respectively.  The  approximate  total  acoustic  field  is  now 


known  from  Eqs.  (10)  and  (11) . 

The  validity  of  the  expansions  in  this  section  and  our 

a 

rationale  for  retaining  second-degree  terms  in  t may  be 
explained  by  examining 
A A . 2 

4>(N,t)  = [u>a^  - (N-I)tt]  + wa^t  + ua^t  . (27) 

Equation  (27)  represents  the  phase  of  the  ray  with  N bottom 

_ A 

reflections  in  terms  of  source  time  t and  follows  from  Eqs.  (17) 

and  (8b) . Since  w is  typically  10^  rad  sec  ^ and  assuming  that 

—2  A 

M may  be  as  large  as  10  , it  can  be  shown  that  if  |t|  < 5 sec, 

A 2 A A 

the  quadratic  term  rna^t  is  small  and  4>(N,t)  may  be  approximated 

A A 

linearly.  For  larger  values  of  |t|  , wa^t^  is  larger  than  M and 

A 

must  be  included  in  the  approximation  of  <j>(Eq.  (27)).  To  determine 

the  time  interval  for  which  Eq.  (27)  is  valid  to  within  a correction 
of  0 (M) , we  derived  the  cubic  term  rna^t  . We  demonstrated  that 
its  magnitude  remains  sufficiently  small  to  assure  the  validity 
of  Eq.  (27)  if  | t | is  < 50  sec  at  Rq  = 5 km  and  for  "large" 
relative  speeds,  although  bigger  time  intervals  may  be  justified 
at  longer  ranges  and  more  moderate  speeds.  Thus,  second-degree 
approximations  improve  the  time  interval  for  which  our  expansions 
remain  valid  while  avoiding  the  inherent  difficulties  of  higher 
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T 


k 


than  second-degree  equations. 

Although  the  previous  paragraphs  are  concerned  with  restrictions 

A 

on  t,  corresponding  limits  on  receiver  time  t may  be  deduced 

A A A A 

from  Eq.  (2a)  (appropriately  replacing  T*(N,t)  with  T (N , t) ) and 
the  fact  that  at  short  ranges  the  ray  with  one  bottom  reflection 
always  exists.  In  particular,  the  difference  between  any 

A 

estimates  of  the  upper  limits  of  t and  t is  the  travel  time  for 
the  N = 1 ray,  which  can  be  shown  to  be  less  than  10  sec  for  the 
ranges  considered  in  this  paper.  Therefore,  the  upper  bound  on  t 

A 

is  slightly  larger  than  for  t,  but  the  discussion  of  restrictions 

A 

on  t is  essentially  valid  for  t,  the  time  variable  of  primary 
interest. 


III.  NUMERICAL  RESULTS  AND  ANALYTICAL  MODEL 

Using  Eqs.  (9b)  - (11)  and  (24) , we  have  obtained  numerical 

results  which  demonstrate  the  effects  on  total-field  amplitude 

and  phase  of  the  variation  of  initial  range  Rq,  receiver  speed 

v^,  source  speed  v^,  and  initial  orientations  represented  by 

the  angles  XQ  and  At  short  ranges,  the  inclination  angle 

of  a ray  at  a bottom  reflection  is  sufficiently  large  that 

bottom  loss  Bq,  or  ratio  of  reflected  to  incident  amplitude, 

17  18 

remains  constant  over  many  ocean  areas.  ' Therefore,  we 


elected  to  use  the  results  of  Ref.  18  which  suggests  Bq  = 0.4  as 


an  appropriate  value.  Also,  we  chose  the  circular  frequency  at 
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} 

remaining  parameter  values  as  follows: 

g^  = .0536  sec  g^  = .0155  sec  \ c^  = 1522  m sec  \ and 
-1 

c^  = 1484  m sec  . We  terminated  the  summations  of  Eqs.  (9b)  - (11) 

at  N = N',  where  N’  + 1 denotes  the  first  ray  to  contribute  less  than 

one  percent  of  the  amplitude  of  the  N = 1 ray.  For  all  parameter 

values  considered  here,  the  summations  could  be  terminated  at  N'  =4. 

Figure  3 illustrates  60-sec  runs  for  R = 2 km  and  10  km  and 

o 

• -1  o 

for  v^  = V2  = 10  kn  (=  5.14  m sec  ) . The  sound  source  and 

receiver  iR,  move  on  the  same  linear  path,  but  in  opposite  directions. 

• -1 

Therefore,  v = 20  kn  (=  10.28  m sec  ),  A - t\  rad,  and  \p  =0  rad. 

o o 

Although  we  show  typical  60-sec  runs,  for  these  particular  receiver 

and  source  speeds  and  directions,  results  based  on  approximations 

such  as  Eq.  (23)  remain  valid  for  approximately  two  minutes  of 

receiver  time  at  Rq  = 2 km,  and  as  discussed  earlier,  for  longer 

periods  of  time  at  R^  = 10  km.  It  can  be  seen  from  Fig.  3 that 

total-field  amplitudes  show  a pseudo-sinusoidal  behavior  in  which 

both  their  extent  and  rate  of  variation  (oscillations  per  sec)  are 

greater  at  initial  range  Rq  = 10  km.  The  total-field  cumulative 

phase  shows  an  overall  parabolic  behavior  in  which  0 < <J>"(t)  <<  1. 

To  demonstrate  this  subtle  curvature,  we  have  added  dashed  straight 

lines  tangent  to  the  actual  cumulative  phase  curves  at  t = 0. 

Fig.  3 illustrates  also  that  this  parabolic  behavior  becomes 

more  pronounced  as  R decreases,  and  at  R =2  km,  there  is  as 

o o 

much  as  an  8-cycle  difference  between  the  dashed  line  and  the  actual 


parabolic  phase  curves.  It  can  also  be  seen  that  the  variation  of 
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phase  over  time,  or  phase  rate,  increases  with  range  for  a given 

relative  velocity  v.  For  example,  at  R = 2 km,  the  variation  is 

o 

approximately  48  cycles  over  a 60-sec  period,  whereas  at  Rq  = 10  km, 
the  variation  is  over  120  cycles.  The  cumulative-phase  values 
shown  at  t = 0 were  chosen  to  be  those  between  zero  and  one  cycle. 

We  add  that  for  runs  in  which  the  receiver-source  distance  decreased 
with  time,  the  cumulative  phase  rate  was  negative,  as  expected, 
and  its  magnitude  varied  with  range  similar  to  that  just  described 
for  positive  phase  rate.  Also,  the  parabolic  nature  of  cumulative 
phase  at  short  ranges  seemingly  contrasts  with  the  essentially 
linear  behavior  of  phase  at  long  ranges,  illustrated  and 
discussed  in  Refs.  3 thru  6.  However,  the  long-range  assumptions 
of  these  studies  permitted  linear  time-dependent  approximations  of 
phase  for  reasonable,  but  suitably  restricted,  time  intervals.  If 
the  restrictions  on  time  in  these  investigations  had  been  relaxed 
sufficiently  to  warrant  second-degree  approximations,  a nonlinear 
pattern  for  cumulative  phase  might  also  be  expected  at  long 
ranges  for  sufficiently  long  sampling  periods.  Finally,  the 
extent  (3-5dB)  and  more  frequent  amplitude  fades,  in  addition  to 
the  well-behaved  amplitude  pattern,  at  the  short  ranges  in  this 
paper  contrasts  sharply  with  the  deeper  (20-25  dB)  and  less 
frequent  fades  of  the  more  noise-like  amplitude  pattern  which 
these  references  have  described  for  longer  ranges. 


The  consistent  and  well-behaved  patterns  exhibited  by 
total-field  amplitude  and  phase  throughout  all  numerical  results 
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suggest  that  it  should  be  possible  to  obtain  approximate 

analytical  expressions  to  explain  the  behavior  of  the  total- 

field,  rather  than  relying  on  numerical  computations  for 

particular  sets  of  parameters.  Therefore,  we  forego  displaying 

additional  numerical  results  and  turn  to  the  development  of  an 

analytical  model.  Using  Eqs.  (13)  and  (25),  it  can  be  shown  that 

-1  -2 

A ( N ) /A ( 1 ) is  0(10  ) for  N = 2,3  and  is  0(10  ) for  N > 3.  With 

these  magnitudes  of  relative  amplitude,  Eq.  (10)  and  basic 
trigonometric  identities  can  be  used  to  write 

CL (t)  = A(l)  + B (t) cos  [n  (t)  ] , (28) 


where 

2 2 1/2 
B(t)  = { [A(2) ] + [A(3)l  + 2A(2) A( 3) cos [uT ( 3 , t) -wT (2 , t) ] } ' , (29a) 

and 

n (t)  = ai[T(2,t)-T(l,t)  ] + tan_1{A(3)sin[a>T(3,t)-u>T(2,t)  ] } (29b) 

x {A  ( 2 ) + A ( 3 ) cos  [idT  (3 , t)  -<dT  ( 2 , t)  ] }_1  , 


where  T(N,t)  is  known  from  Eq.  (23)  and  A (N)  from  Eq.  (26a). 

To  approximate  $(t),  we  proceed  in  a manner  somewhat  similar 
to  derivations  contained  in  Ref.  19.  Specifically,  if  ACL(t)  and 
A4>(t)  represent  the  total-field  deviations  due  to  motion,  then 


o a(t)  = a(t)  - <a. 


and 


M(t)  = 4>(t)  - $> 


(30a) 


(30b) 


where  and  4>q  are  total-field  amplitude  and  phase,  respectively, 

without  receiver  and  source  motion.  Using  Eqs.  (30) , Eq.  (9b) 
can  be  written  as  the  imaginary  part  of 

(&0+ACl)exp{i  [o)t-4>o-A4>]  } = Z A (N)  exp{  i [ut-$  (N,  t)  ] } . (31) 

N 

If  both  sides  of  Eq.  (31)  are  multiplied  by  exp  {-iltot  - 4>q]  } , 
and  the  terms  of  the  summation  on  the  right  are  combined  using 
basic  trigonometric  identities,  it  can  be  shown  that 

*(t)  = <i)T  ( 1 , t)  + [A(l)]_16(t)sin[n(t)]  , (32) 

-2 

where  we  drop  terms  of  size  less  than  10  , AG./CL  <<  1,  and 

o 

A(N)/A(1)  diminishes  with  N,  as  discussed  previously. 

Comparisons  of  the  numerical  results  obtained  from  Eqs.  (9b)  - 
(11)  and  (26) , and  the  analytical  approximations  represented  by 
Eqs.  (28)  and  (32) , show  that  Eq.  (28)  consistently  approximates 
CL(t)  to  within  0.5  dB,  and  Eq.  (32)  approximates  4>  ( t ) to  within 
0.05  cycles.  In  all  cases,  receiver  time  was  restricted  according 
to  initial  range  and  relative  speed,  as  previously  discussed.  For 
plotting  purposes  now  and  to  simplify  subsequent  analyses,  the 
right  side  of  Eq.  (32)  may  be  approximated  by 

<t>  (N,  t ) = wlt+fa^-t)  d+a2)  1 + a3  (a^^-t)  2 (l+a2)  3]-(N-1)tt  (33) 

evaluated  at  N=l,  where  Eq.  (33)  follows  from  substituting  Eq.  (23) 
into  Eq.  (26b).  Figure  4 compares  numerical  and  analytical 
values  for  a magnified  portion  of  the  numerical  results  in  Fig.  3. 


The  results  shown  are  representative  of  the  entire  60-sec  run.  The 


solid  curves  are  numerical  results,  the  dashed  amplitude  curves 
are  graphs  of  Eq.  (28)  , and  the  dashed  cumulative  phase  curves 
depict  Eq.  (33)  evaluated  at  N=l.  The  approximate  phase  values 
corresponding  to  the  complete  Eq.  (32) , which  contains  multipath 
contributions,  generally  fall  between  the  solid  and  dashed  phase 
curves  of  Fig.  4.  Thus,  cumulative  phase  may  be  conveniently 
described  as  the  sum  of  the  dominant  phase  of  the  N=1  ray, 
approximated  by  Eq.  (33) , and  a multipath  correction  due  to  the 
combined  contributions  of  the  rays  corresponding  to  N > 2.  This 
correction,  approximated  by  the  second  term  of  the  right  side  of 
Eq.  (32),  causes  the  total-field  phase  to  oscillate  subtly  about 
the  phase  of  the  N=1  ray.  This  tendency,  confirmed  analytically 
here,  was  observed  also  in  the  computer  simulations  in  Ref.  3. 

Finally,  in  contrast  to  the  excellent  single-ray  approximation  to 
cumulative  phase,  Eq.  (28)  shows  that  the  approximation  to  total-field 
amplitude  relies  on  the  effects  of  multipath. 


IV.  TOTAL- FI ELD  STUDY 

We  turn  next  to  an  examination  of  the  total  field  using  the 
analytical  model  developed  in  Sec.  III.  This  section  consists  of 
quantitative  descriptions  of  effects  of  the  variation  of  initial 
range  and  changes  in  the  parameters  of  motion  on  total-field 
amplitude  and  phase.  For  brevity,  we  have  omitted  other  results 
which  could  be  obtained  from  Eqs.  (28)  and  (32).  As  examples. 


it  is  possible  to  describe  the  relationship  between  amplitude 
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oscillations 


and  cumulative  phase. 


as  well  as  to  show  that  the 


extent  of  amplitude  fades  increases  with  range  but  is  independent 
of  a nonzero  relative  speed. 

For  purposes  of  discussion,  we  group  the  effects  of  receiver 
and  source  motion  into  two  categories.  The  first  category  consists 
of  effects  attributable  to  the  relative  motion  between  the 
receiver  and  source,  represented  by  terms  containing  the 
parameter  M (or  equivalently,  v) . If  the  receiver  is  stationary, 
only  this  category  exists  and  the  relative-motion  parameters  are 
determined  completely  from  those  of  source  motion.  However,  if 
the  receiver  is  in  motion,  there  is  a second  category  of  effects 
due  to  an  apparent  current,  represented  by  terms  containing  the 
parameter  (or  v^  = v^) . Furthermore,  since  a moving  receiver 
also  affects  relative  motion  (see  Eqs.  (3c)  and  (3d)),  receiver- 
motion  contributions  to  the  total  field  are  included  in  both  of 
the  groups  just  described.  Only  the  second  category  (current 
effects)  is  considered  now,  while  that  dealing  with  relative 
motion  will  be  addressed  subsequently. 

As  demonstrated  earlier,  Eq.  (33)  evaluated  at  N=1  closely 
approximates  total-field  cumulative  phase.  If  Eqs.  (18),  (19a), 

and  (20a)  are  substituted  into  Eq.  (33)  and  if  we  ignore  the 
small  terms  for  which  the  sum  of  the  exponents  of  M and  M 


exceeds  two,  then 


4>  = 0){T  (N)  [ 1-M  cos0,  cosX  - MG  c R cos  \1)  ) 
o c bo  o 2 b o ro 


+ ‘WbC  + HG2°bRolcOS  *o]t 

+ M2H  (c.  R_1 ) 2 1 2 } | 

^ O I M l 


where  G , G^,  and  H are  given  by  Eqs.  (1%),  (lgc),  and  ( 2 0 ) , 

respectively.  We  have  demonstrated  numerically  that  Eq.  (34) , 

although  simpler  than  Eq.  (33) , is  an  excellent  approximation  for 

total-field  phase.  We  have  also  shown  that  G. c.  R 1 is  0(1),  and 

I b o 

G.c.  R 1 is  0(10  ^)  for  R < 3 km  and  0(1)  for  3 < R <10  km. 

Z D O O ~ ~ O — 

-12.  -2  -1 
Also,  H ( c^R^  ) increases  from  approximately  10  sec  to 

10  ^sec  ^ as  Rq  decreases  from  10  km  to  2 km.  Thus,  Eq.  (34) 

shows  that  current  contributions,  or  second-category  effects. 


are  greater  than  M for  both  the  time- independent  and  time- 

dependent  portions  of  total-field  phase.  Also,  in  certain  cases, 

these  effects  dominate  those  of  relative  motion.  In  particular, 

• 3 -1 

for  a)  = 10  rad  sec  , the  time-independent  current  effect 
represented  by  cos0^ocosXo  ranges  from  approximately 

one  to  ten  radians,  depending  on  range,  for  cos  Xq  not  close 
to  zero.  Suppose,  also,  that  the  component  of  relative  velocity 
vcos  is  close  to  zero  ; i.e.,  either  the  path  of  source  motion 


is  essentially  perpendicular  to  the  receiver-source  line  or 
relative  speed  is  negligible.  Then  the  time-independent  current 
contribution  dominates  that  of  relative  motion,  represented  by 
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the  term  wT  ( 1) MG_c.R  cos  U/  . If  cos  \ is  close  to  zero, 

o 2 b o o o 

however,  the  current  is  approximately  perpendicular  to  ray-path 

planes,  and  its  effects  are  obviously  minimal.  Similarly, 

Eq.  (34)  shows  that  for  most  receiver  and  source  directions,  the 

time-dependent  current  contribution,  represented  by  the  term 

oiM^MG^c^R^^t , is  small  in  comparison  to  the  time-dependent 

effects  of  relative  motion.  However,  it  is  easy  to  show  that  the 

variation  of  total-field  phase  due  to  current  can  be  as  much  as 
o 

200  over  a two-minute  time  period  for  receiver  and  source  speeds 
of  approximately  10  kn  and  source  frequency  of  10^  rad  sec  ^ . 

Thus,  if  either  the  path  of  relative  motion  is  virtually 
perpendicular  to  receiver-source  direction  or  relative  speed 
is  negligible,  the  time  dependence  of  total-field  phase  can  be 
attributed,  within  our  approximations,  to  current  effects. 

In  order  to  simplify  our  analysis  of  effects  described  earlier 
as  belonging  to  the  first  category  ( relative  motion  effects), 
we  introduce  a time-independent  approximation  to  phase  rate  4> ' . 
Using  Eqs . (19a)  and  (34),  we  may  approximate 

with  the  N=1  ray  as  follows: 


= coG_c  R M cos  \p 
2 b o o 


= w(c°s9boCscemo[sinemo-sineboJ [2(g11+g21)l 

X ♦ V,4(3i1+g21)cbI)) 


X R COS  ill 

O O > 
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where,  based  on  previous  discussions  of  coefficient  magnitudes, 

2 

we  have  ignored  the  McMt  and  M t terms.  We  remark  that  the 
approximation  for  4> ' could  be  improved,  for  example,  by  using  an 
average  over  time,  or  possibly  other  procedures.  For  simplicity, 
we  proceed  using  Eq.  (35) , which  can  be  used  to  demonstrate 
effects  of  relative  motion  on  total-field  amplitude  and  phase 
over  limited  time  intervals.  We  note  also  that  since  4>(t)  generally 
is  not  proportional  to  time,  the  concept  of  frequency  is  not  well 
defined;  however,  our  time-independent  approximation  of  4> ' may  be 

interpreted  as  an  estimate  of  Doppler-shif t frequency  with  units  of 

. -1  6,20 

rad  sec 

In  Fig.  5,  the  straight  lines  (solid  for  Rq  = 2 km  and  10  km, 
and  dashed  for  5 km)  follow  from  Eq.  (35) . They  show  effects  for 
different  ranges  of  changes  in  the  component  of  relative  speed  in 
the  initial  receiver-source  direction,  v|cos\J>o|,  on  a normalized 
phase  rate  w ^3>'.  Figure  5 illustrates  that  the  effect  on  phase 
rate  of  differences  in  relative  speed  is  reduced  as  initial 
range  Rq  decreases,  and  that  for  a given  component  of  relative 
speed,  differences  in  range  correspond  to  smaller  differences  in 
phase  rate  at  longer  ranges.  This  suggests  that  as  range  increases 
beyond  the  approximate  10  km  limit  in  this  paper,  effects  of 


range  differences  are  reduced,  and  differences  in  the  relative- 


speed  component  become  the  dominant  factor  in  determining  phase 
rate.  Finally,  as  expected,  phase  rate  is  positive  (negative) 
according  to  receiver-source  distance  increasing  (decreasing) 
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with  time,  or  equivalently,  vcos^q  being  positive  (negative) . 

Moreover,  if  vcosil*  = 0,  then  the  effects  of  relative  motion  on 

o 

phase  rate  are  negligible,  although  current  effects  may  cause 

phase  to  exhibit  some  time  dependence,  as  discussed  earlier.  To 

demonstrate  the  usefulness  of  Eq.  (35)  as  an  estimate  of  phase 

rate  for  a specified  time  interval,  if  oi  = 2ir(350)rad  sec  ^ and 
• -1 

vcos\Jjq  = 2 kn  (=  10.28  m sec  ),  then  from  Fig.  5 at  Rq  = 2 km, 

• 3 

♦ ' = (2.042)a)  • 10  = 0.72  Hz.  This  value  compares  favorably 

with  the  phase  rate  of  0.80  Hz,  shown  in  the  numerical  results  of 
Fig.  3 as  a total  change  in  phase  of  approximately  48  cycles  over 
a 60-sec  time  period. 

Equations  (28)  and  (29)  can  be  used  to  describe  how  the  time 
between  successive  amplitude  peaks,  which  we  denote  by  y,  varies 
with  relative  motion  at  different  ranges.  In  particular,  Eq.  (28) 
shows  that  CL  ( t)  reaches  a relative  maximum  with  B ( t) cos  [n (t) ] . 
Also,  Eq.  (29b)  indicates  n(t)  = <j>(2,t)  - <j>(l,t),  since  A(3)/A(l) 
is  0(10  1)  and  therefore,  cos[n(t)]  peaks  when 
cos[(|>(2,t)  ~ <M1,  t)  ] = 1.0.  Approximating  B(t)  with  A(2),  we 
conclude  that  8 (t) cos [n (t) ] peaks  at  approximately  the  same  time 
as  A(2)  cos  [<f>  (2 , t)  - 4>(l,t)].  Also,  if  4>'(2,t)  is  approximated 
linearly  in  the  same  manner  as  <j>'(l,t)  in  Eq.  (35),  we  obtain 

n (t)  = <M2,t)-<Ml,t)  = a)  [G_  (2) -G  (1)  ] R v(cosiJ>  ) t+K  , (36) 

z 2 O O 

where  K is  a constant.  From  Eq.  (36) , it  follows  that 

(37) 
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where  f is  frequency  in  Hz.  Equation  (37)  provides  a simple 

expression  with  which  to  analyze  the  effects  on  p of  changes 

in  relative-speed  component  and  range.  The  hyperbolas  in  Fig.  5 

illustrate  such  effects  for  a normalized  peak-to-peak  time 

pf  at  Rq  = 2 km  and  10  km  (solid  curves)  and  at  an  intermediate 

range  of  5 km  (dashed  curve).  In  contrast  to  phase  rate,  which 

essentially  increases  or  decreases  proportionally  to  vcos^, 

differences  in  peak-to-peak  time  pf  corresponding  to  differences 

in  relative-speed  component  are  much  greater  at  lower  than  at 

higher  speeds.  Also,  for  a particular  value  of  v]cos\|>o|,  a 

change  in  range  corresponds  to  a greater  change  in  pf  at  lower 

speeds;  and  similar  to  the  variation  in  phase  rate,  differences 

in  range  correspond  to  smaller  differences  in  pf  at  longer  ranges. 

As  expected,  if  vjcos^^l  is  small,  pf  becomes  large,  consistent 

with  amplitude  becoming  constant  over  time  as  relative  motion 

becomes  negligible.  In  addition,  the  results  using  Eq.  (37)  and 

shown  in  Fig.  5 agree  favorably  with  numerical  results.  For 

• -1 

example,  if  v = 20  kn  (=  10.28  m sec  ),  Rq  = 2 km,  f = 350  Hz, 
and  cosi(io  = 1,  then  from  Fig.  5,  p = 2.86  sec.  For  the  60-sec 
run  in  Fig.  3 there  are  approximately  23  fades  in  60  sec  for  an 
average  time  between  amplitude  peaks  of  2.61  sec.  Figure  5 
suggests  also  that  for  a given  relative-speed  component,  the 
time  between  amplitude  peaks  achieves  a relative  minimum  in 
the  interval  2 km  < Rq<  5 km.  This  can  be  confirmed  analytically. 
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a function  of  initial  range  Rq,  has  a relative  maximum  in  the 
same  interval. 

As  noted  earlier,  acoustical  effects  of  a moving  receiver 

can  be  classified  as  those  due  to  a uniform  current,  already 

discussed,  and  those  resulting  from  changes  to  relative 

receiver-source  motion.  In  order  to  consider  the  latter,  it  is 

necessary  to  examine  first  the  relative-motion  parameters  as 

functions  of  receiver-motion.  In  the  trivial  case  of  a fixed 

source  (v2  = 0)  , receiver  motion  completely  determines  the 

relative  motion  and  the  effects  are  as  described  in  the  preceding 

paragraphs.  If  / 0,  the  effects  of  a moving  receiver  on 

relative  motion  can  be  portrayed  as  in  Figs.  6(a)  and  6(b), 

using  Eqs.  (3c)  and  (3d),  respectively.  Figure  6(a)  shows  the 

effect  of  changes  in  a dimensionless  receiver  speed  v^/v^  on  a 

dimensionless  relative  speed  v/v^,  while  Fig.  6(b)  shows  the  effects 

of  changes  in  v^/v2  on  cos^q.  For  each  figure,  the  directions  of 

source  and  receiver  motion  at  time  t=0  are  represented  in  the 

term  cos(Xq-yo),  whose  values  generate  a family  of  curves 

describing  the  possible  relationships  illustrated.  Only  the  solid 

portions  of  the  parabolas  in  Fig.  6(a)  are  valid,  since  speeds 

are  nonnegative.  It  can  be  seen  that  for  cos (A  -v  ) < 0,  v/v 

o o — 2 

increases  with  v^/v^  an<^  s9nlc°s\f>0]  is  independent  of  v^/v^. 

However,  if  cos(Ao~yo)  > 0,  v/v 2 decreases  before  increasing 
with  v-^/V2  an<^  sgn  [cosi^l  depends  on  v^/v2 . In  addition,  for 
Fig.  6(b),  if  cos(xo-yq)  / + 1,  there  are  two  possible  curves 
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corresponding  to  cos^  > 0 and  cos^q  < 0,  respectively;  i.e., 

the  relative  motion  of  the  source  at  t=0  may  be  either  away  from 

or  toward  the  receiver.  Finally,  there  are  two  possibilities 

each  for  cos  (A  -y  ) = + 1.  If  cos  (A  -y  )=-l,  then  cosy  = +1(-1) 
o o — o o o 

for  all  vj/v2'  if  cosyQ  = +1(-1);  i.e.,  receiver-source  distance 

is  increasing  (decreasing)  with  time.  However,  if  cos (A  -y  ) = 1 

o o 

and  cosyQ  = 1,  then  cosi|>q  jumps  from  +1  to  -1  at  v^/v  = 1*  if 
cosYo  = -1,  cosi/>q  jumps  from  -1  to  +1  at  v^/Vj  = T^us'  v^/v^ 

determines  the  sign  of  the  rate  of  change  of  receiver-source 

distance  as  was  the  case  for  cos (A  -y  ) > 0. 

o o 

Equations  (3c)  and  (3d)  can  be  substituted  into  Eqs.  (35)  and 
(37)  in  order  to  examine  specifically  the  effects  of  receiver 
motion  on  total-field  amplitude  and  phase  for  a particular  source 


speed  artd  range.  As  an  illustration,  we  show  such  effects  for 
• -1 

v = 20  kn  (=  10.28  m sec  ) and  R = 2 km  in  Figs.  7 and  8. 

« o 

It  can  be  seen  from  Fig.  7 that  for  nonpositive  values  of 

cos (A  -y  ),  since  the  time  rate  of  change  of  receiver-source 
o o 

distance  does  not  change  sign  and  either  increases  or  decreases 

with  v^,  the  correspondence  between  u)  and  receiver  speed  is 

straightforward  and  virtually  linear.  However,  if  cos (A  -y  ) > 0, 

o o 

differences  in  w due  to  differences  in  receiver  speed  depend 

greatly  on  the  magnitude  of  v^,  including  a possible  difference 

in  sign.  For  example,  if  cos (A  -y  ) = 1/2  and  the  source  is 

o o 

moving  away  from  the  receiver,  there  is  an  approximate  difference 


31 


of  1.0  in  10  a)  4>  ' corresponding  to  the  two  values  v^/v^  = ^ 

and  0.5.  However,  if  the  relative  speed  is  between  0.5  and  1.5, 

the  corresponding  difference  in  lo\)  ' is  small;  if  v^/v2 

exceeds  2.0,  lo\i  ^<t> ' is  negative.  This  pattern  of  behavior  of 

a)  14>*  is  a direct  result  of  the  earlier  discussions  regarding 

the  effects  of  changes  in  receiver  speed  v,  on  both  v and  cosiJj 

1 o 

when  cos  (Xq-yo)  > 0»  which  can  be  seen  in  Figs.  6(a)  and  6(b) 
for  the  same  values  of  V^/V2* 

Trends  for  yf  are  shown  in  Fig.  8,  using  the  same  values  of 

Rq  and  source  speed  v2  as  in  Fig.  7.  As  discussed  previously, 

if  cos  (Xq  - yq)  < 0,  the  time  rate  of  change  of  receiver-source 

distance  does  not  change  sign,  and  its  magnitude  increases  with 

v^.  This  is  equivalent  to  |vcos^q|  being  bounded  away  from  zero 

and  increasing  with  v^.  The  result  is  that  yf  decreases,  but  at 

a rate  which  depends  on  the  receiver  and  source  directions,  or 

on  cos  (X  - y ).  However,  if  cos(X  - Y ) >0,  the  variation 
o o o o 

of  yf  with  receiver  speed  v^  is  less  predictable,  as  was  the  case 
for  the  variation  of  w with  v^.  Also,  the  tendency  of  yf  to 

increase  without  bound  corresponds  to  total-field  amplitude 
becoming  constant.  This  geometric  consequence  can  be  attributed 
to  that  particular  ratio  of  receiver  and  source  speeds  which,  for 
a given  combination  of  receiver  and  source  directions,  causes  the 
relative-speed  component  to  become  zero  and  results  in  a negligible 
time  rate  of  change  of  the  horizontal  distance  traveled  by  the 

k ■ =r  - 
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sound  rays.  Also  from  Fig.  8,  if  the  receiver  is  stationary 

-3  • -1 

(v^/v^  = 0),then  10  pf  = 1.  Since  v ^ = 20  kn  (=  10.28  m sec  ) 

and  Rq  = 2 km,  if  f = 350  Hz,  then  p = 2.87  sec.  This  agrees 

closely  with  the  numerical  results  of  Fig.  3 (discussed  earlier) 

for  which  both  receiver  and  source  were  moving,  but  with  the 

same  relative  speed  of  20  kn.  It  can  be  seen  from  Fig.  8, 

however,  that  for  > 0,  the  effects  of  receiver  motion  are 

distinctly  different,  depending  on  the  receiver  and  source 

directions.  We  note  that  graphs  similar  to  those  in  Figs.  7 

and  8 could  be  sketched  for  other  parameter  choices,  although 

the  trends  would  remain  the  same.  Also,  for  smaller  values  of 

source  speed,  the  intercepts  on  the  vertical  axes  decrease  for 

o>  and  increase  for  pf,  and  conversely  for  larger  values. 

Such  figures  illustrate  vividly  the  significant  effects  on  total 

field  which  can  be  caused  by  differences  in  receiver  speed, 

depending  on  receiver  and  source  directions,  source  speed,  and 

range. 

V.  SUMMARY 

This  paper  presents  an  analysis  of  the  effects  of  source  and 
receiver  motion  on  a CW  signal  transmitted  in  a deep-ocean  channel 
over  short  ranges.  The  source  and  receiver  each  follow  constant- 
velocity  horizontal  paths  on  the  ocean  surface.  A bilinear  sound- 


speed  profile  and  horizontal,  smooth  boundaries  are  assumed.  The 
moving  source,  moving  receiver  problem  is  transformed  into  a 
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geometrically  and  acoustically  equivalent  fixed-receiver  model. 

Under  certain  time  restrictions,  general  equations  are  derived 
for  the  launch  angles,  travel  times,  and  spreading  losses  of 
SRBR  rays,  which  lead  to  time-dependent  total-field  phase  and 
amplitude  expressions.  Because  of  the  short  ranges,  linear 
approximations  in  previous  studies  are  replaced  by  second- 
order  approximations  in  order  to  maintain  mathematically  valid 
expressions  over  time  intervals  of  sufficient  length. 

Total-field  amplitude  and  cumulative  phase  are  calculated 
as  functions  of  time  for  various  source/receiver  speeds, 
orientations,  and  ranges.  It  is  shown  that  amplitude  displays 
a rapidly  oscillating  pseudo-sinusoidal  behavior  and  that  cumulative 
phase  exhibits  an  overall  parabolic  pattern,  with  a subtle 
oscillatory  behavior  on  a smaller  scale.  These  results  contrast 
with  the  more  noise-like,  but  slower  oscillating,  amplitudes 
and  the  essentially  linear  behavior  of  cumulative  phase  over 
time  demonstrated  in  previous  investigations  for  comparable  time 
intervals  but  at  longer  ranges. 

Analytical  expressions  are  derived  which  closely  approximate 

» 

numerical  results.  These  expressions  provide  the  basis  for 
studying  total-field  amplitude  and  phase  in  terms  of  arbitrary 
motion  parameters  and  initial  range.  It  is  shown  that  phase  rate 
increases  strictly  with  both  range  and  the  component  of  relative 
speed  in  the  receiver-source  direction,  with  the  effects  of 
changes  in  range  diminishing  at  larger  ranges.  Peak-to-peak 
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amplitude  time  is  shown  to  decrease  with  an  increase  in  relative 
speed,  with  the  effects  of  changes  in  speed  diminishing  at 
higher  speeds  and  the  effects  of  differences  in  range  diminishing 
with  higher  ranges. 

The  effects  of  receiver  motion  on  total  field  are  illustrated 
and  discussed.  It  is  demonstrated  that  a moving  receiver 
contributes  to  the  total  field  by  introducing  a current  effect 
and  also  by  changing  the  relative  motion.  In  general,  the 
primary  current  contribution  is  a constant  shift  in  total-field 
phase,  with  the  dependency  of  phase  on  time  being  dominated  by 
relative  motion.  However,  for  certain  receiver-source  directions, 
the  time-dependent  current  effect  can  cause  significant  variations 
in  total-field  phase  and  the  time-dependence  of  phase  may 
essentially  be  attributed  to  only  current  effects.  Significant 
receiver-motion  effects  on  total  field,  through  changes  in  relative 
motion,  are  demonstrated,  to  include  the  relative  importance  of 
receiver-source  directions,  as  well  as  the  ratio  of  receiver  speed 
to  source  speed. 


35 


P 

REFERENCES 

R.B.  Lindsay,  Mechanical  Radiation  (McGraw-Hill,  New  York 
1960),  pp.  311-321. 


o 

P.M.  Morse  and  K.U.  Ingard,  Theoretical  Acoustics  (McGraw-Hill, 
New  York,  1968),  pp.  699-726. 


^R.P.  Flanagan,  N.L.  Weinberg,  and  J.G.  Clark,  J.  Acoust.  Soc. 
Am.  56,  1673-1680  (1974). 


4 

J.G.  Clark,  R.P.  Flanagan,  and  N.L.  Weinberg,  J.  Acoust.  Soc. 
Am.  60,  1274-1284  (1976). 


^G.M.  Jacyna,  M.J.  Jacobson,  and  J.G.  Clark,  J.  Acoust.  Soc.  Am. 
60,  815-824  (1976). 


^G.M.  Jacyna  and  M.J.  Jacobson,  J.  Acoust.  Soc.  Am.  _61,  1153-1162 
(1977) . 

7 

J.A.  Neubert,  J.  Acoust.  Soc.  Am.  62,  1404-1411  (1977). 


^R.E.  Wilcox,  J.  Acoust.  Soc.  Am.  63,  870-875  (1978). 


Q 

G.M.  Jacyna  and  M.J.  Jacobson,  J.  Acoust.  Soc.  Am.  6J3  , 1353-1364 
(1978)  . 


■^R.J.  Urick,  Principles  of  Underwater  Sound  (McGraw-Hill),  New  York, 
1975,  2d  edition),  pp.  53-54. 


^C.B.  Officer,  Introduction  to  the  Theory  of  Sound  Transmission 
(McGraw-Hill,  New  York,  1958),  pp.  228-248. 


12R.N.  Baer  and  M.J.  Jacobson,  J.  Acoust.  Soc.  Am.  54,  90-91  (1973). 


13 

J.A.  Widtfeldt  and  M.J.  Jacobson,  J.  Acoust.  Soc.  Am.  59 , 
852-860  (1976). 


14 

K.G.  Hamilton,  W.L.  Siegmann,  and  M.J.  Jacobson,  J.  Acoust. 

Soc.  Am.  62,  53-62  (1977). 

"'■^L.A.  Stallworth  and  M.J.  Jacobson,  J.  Acoust.  Soc.  Am.  S2,  344-355 
(1972)  . 

18B.K.  Newhall,  M.J.  Jacobson,  and  W.L.  Siegmann,  J.  Acoust.  Soc. 

Am.  62,  1165-1175  (1977). 

17Ref . 10,  pp.  126-132. 

18K.V.  Mackenzie,  J.  Acoust.  Soc.  Am.  32^,  221-231  (1960)  . 

19 

W.L.  Siegmann,  M.J.  Jacobson,  K.G.  Hamilton,  "Effects  of  Tidally- 
varying  sound-speeds  in  acoustical  propagation  over  a sloping 
ocean  bottom,"  Rensselaer  Polytechnic  Institute,  Math.  Rep. 

No.  123,  1 February  1979  (unpublished). 


2°Ref . 2,  pp.  717-732. 


FIGURE  LEGENDS 


FIG.  1 


FIG.  2 


FIG.  3 


FIG.  4 


FIG.  5 


FIG.  6 


FIG.  7 


FIG.  8 


(a)  Bilinear  sound-speed  channel  and  (b)  geometry 
for  constant-velocity  receiver  and  source  motions. 

Moving  reference  frame  geometry  for  receiver  and 
source  motions  shown  in  Fig.  1. 

Cumulative  phase  and  amplitude  vs  time  over  a 60-sec 

interval  for  R = 2 km  and  10  km  with  v,  = = 10  kn, 

o 12 

iji  = 0 . 

o 

Approximate  (dashed)  and  numerical  (solid)  cumulative 
phase  and  amplitude  for  portion  of  run  in  Fig.  3. 

Parameter  values  as  in  Fig.  3. 

Normalized  phase  rate  (lO^u)  "*‘4>')  and  normalized  peak- 
to-peak  amplitude  time  (10  ^pf)  vs  relative-speed 
component  (v|cos\Jjo|)  at  different  initial  ranges. 

(a)  Dimensionless  relative  speed  (v/v^)  and  (b)  cos^  vs 

dimensionless  receiver  speed  (v^/v^)  f°r  different  receiver 

and  source  directions  (cos(X  - y )),  where  v_  ^ 0. 

o o 2 

Normalized  phase  rate  (u>  *4>')  vs  dimensionless  receiver 

speed  (v^/v^)  f°r  RQ  = 2 km,  V2  = kn,  an<^  different 

receiver  and  source  directions  (cos(X  — y )). 

o o 

Normalized  peak-to-peak  amplitude  time  (yf)  vs  dimensionless 

9 

receiver  speed  (v^/v 2)  for  Rq  = 2 km,  v^  - 10  kn,  and 
different  receiver  and  source  directions  (cos(X  - y )). 


CUMULATIVE  PHASE  (CYCLES)  , R = 2 KM 


CUMULATIVE  PHASE  (CYCLES)  , Rfl  = 10  KM 


= 2 KM 
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